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Abstract 

In this paper, we estimate the gamma-ray and neutrino fluxes coming from dark matter annihilation 
in a Milky Way framework provided by a recent N-BODY HORIZON simulation. We first study the 
characteristics of the simulation and highlight the mass distribution within the galactic halo. The 
general dark matter density has a typical r~ 3 power law for large radii, but the inner behaviour is poorly 
constrained below the resolution of the simulation (~ 200 pc). We identify clumps and subclumps and 
analyze their distribution, as well as their internal structure. Inside the clumps, the power law is 
rather universal, r -25 in the outer part with again strong uncertainties for smaller radii, especially for 
light clumps. We show a full-sky map of the astrophysical contribution to the gamma-ray or neutrino 
fluxes in this N-body framework. Using quite model independent and general assumptions for the high 
energy physics part, we evaluate the possible absolute fluxes and show some benchmark regions for 
the experiments GLAST, EGRET, and a km3 size extension of ANTARES like the KM3NeT project. 
While individual clumps seem to be beyond detection reach, the galactic center region is promising 
and GLAST could be sensitive to the geometry and the structure of its dark matter distribution. 
The detection by a km3 version of ANTARES is, however, more challenging due to a higher energy 
threshold. We also point out that the lack of resolution leaves the inner structure of subhalos poorly 
constrained. Using the same clump spectrum and mass fraction, a clump luminosity boost of order ten 
can be achieved with a steeper profile in the inner part of the sub-halos. 



1 Introduction 

An ever increasing number of observational and theoretical results strongly suggest, or even require the 
existence of dark matter, whose enigma becomes thus crucial for the understanding of our universe. 
Let us mention amongst others the WMAP results on CMB p], the rotation curves of disk galaxies 
[2], the formation of large scale structures [3], the bullet cluster observation (4], merger modeling and 
lensing results e.g. [§]. Finally, the possibility of numerous extensions of high energy physics beyond 
the standard model (BSM) to provide new weakly interacting massive particle (WIMP) candidates for 
dark matter makes the hypothesis very appealing. 
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Nevertheless, the nature, the identification and the distribution of the dark matter are still open 
questions, intimately linked with the proof of its existence. Present and near future instrumental 
projects could bring welcome input to these questions. Namely, both astroparticle physics and collider 
searches will reach higher sensitivities with experiments like LHC (accelerator); EDELWEISS II, su- 
perCDMS and ZEPLIN (direct detection); ANTARES and ICECUBE (neutrino telescopes); GLAST 
(gamma space telescope) and PAMELA (charged particle search satellite). Amongst the different 
possibilities, indirect detection is particularly promising. Indeed, relic dark matter particles can accu- 
mulate in cosmic storage rings and annihilate. The decay of their annihilation products will give rise 
to secondary particle fluxes (7, f, e + ,p), which could be detected by dedicated experiments indirectly 
indicating the presence of dark matter. 

In this paper, we will focus on indirect detection of dark matter through gamma rays and neutrinos. 
Galaxies are thought to be interesting sources for this kind of detection, seen the amount of dark matter 
they are believed to harbour. As we will discuss later, the detectability of such gamma rays or neutrinos 
depends strongly on both the astrophysical assumptions on the dark matter distribution in the halo 
and on the assumed high energy physics BSM scenario. Some studies concerning different particle 
physics models can be found in the literature (see (6] [7] for reviews). The popular BSM dark matter 
scenarios are typically supersymmetric models, models with extra dimensions, light dark matter, little 
Higgs model, inert doublet model ... or any extensions providing WIMP. The Milky Way astrophysical 
framework is commonly simplified with assumptions of spherical symmetry, now known to be incorrect 
[3 El [IHl E], and typical smooth dark matter density functions extracted from N-body simulations 
[121 [13] [14] . Few recent works [HI HU [17] and especially [18] with an impressive resolution treat in 
detail the astrophysical aspects of the gamma ray fluxes coming from dark matter annihilation in 
realistic simulation frameworks. Other works e.g [19]. [20] |2"T] [22] [23] consider also more sophisticated 
parametrization inspired by extrapolations of simulation results. 

Typical simulation canvas consist of 1-100 millions of particles with mass around 10 5 — 10 6 solar 
masses. The results reproduce well the large-scale structure formation and have now shown that 
virialized systems are still left with surviving subhalos, also called clumps . The results are more and 
more promising with computing upgrades and resolution improvements. Nevertheless, some questions 
are still open with regard to observations. For instance, the radial density profiles predicted for the 
innermost region of galactic halos are quite cuspy, whereas observations suggest flat cores (see [2j for 
a review). Furthermore, simulations predict more numerous galactic satellites than observed for the 
Milky Way. Even if N-body calculations may generate too concentrated objects, the simulated haloes 
are the only realistic or advanced dark matter distribution framework. Specifically, the estimation of 
dark matter detectability in our neighborhood depends on both the dark matter distribution in the 
Milky Way - especially in the innermost region - and on the number, the size and the concentration of 
the clumps. Depending on the assumptions or results on these key points, different results have been 
proposed in previous works pTl[2l][2^[2^[271[2^[2^[29l[23l[30]. 

The present article is devoted mainly to the astrophysical contribution concerning dark matter 
gamma and neutrino indirect detection. We calculated the possible gamma and neutrino fluxes sky 
map in a N-body simulation framework provided by a HORIZON project simulation [31]. The paper 
is organized as follows: section 2 gives the analysis of the numerical simulation and highlights the 
resulting dark matter distribution. In section 3, the gamma ray and neutrino flux calculation is shortly 
reviewed and a comparison of our estimates with regard to GLAST and ANTARES reach is presented. 
Conclusion and perspectives are given in section 4. 

2 Simulation characteristics 
2.1 General features 

The data used for this paper were provided by the Horizon collaboration. The simulation was performed 
using the Adaptive Mesh Refinement code RAMSES [32]. The initial conditions are set by the WMAP3 
results (Q m = 0.24, fi A = 0.76, Q b = 0.042, n = 0.958, H = 73, cr 8 = 0.77) and the effective number 
of particles is N p — 1024 3 in a box of size L = 20/i _1 Mpc. At z — we selected a Milky- Way sized 
halo. Using the so-called "zoom" technique, we re-defined the grid outside a sphere of diameter 5h , 
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Figure 1: View of the galactic dark matter halo in the HORIZON simulation (projected along the z axis). The 
box size corresponds to 750 kpc. 



using high-mass particles to sample the large scale tidal field and smaller ones for the selected halo 
region. In the high-resolution region, we increase the resolution of the grid on a cell-by-cell basis, with 
a maximum of 7 additional levels of refinement, corresponding to a maximum linear resolution of about 
200 pc. Cells are refined if the local number of dark matter particles exceeds 10. Our smallest particle 
mass is M p — 7.46 10 5 in solar mass (M sun ) units. Our simulated galactic environment is depicted in 



The size of a galactic halo is characterized by its virial radius, r V ir, often defined as the size of 
the sphere centered on the galaxy center with an average density equal to 200 times the cosmological 
matter density. In our simulation, the virial radius is equal to 253 kpc, corresponding to an enclosed 
mass of 6.05 x lO n M or 8.1 x 10 5 particles. 

Dark matter halos are often parameterized by spherically symmetric profiles of the form 



where po is the local density in the solar neighborhood, r$ = 8 kpc is the distance from the sun to the 
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Figure 2: Density at all the positions of the particles in the simulation as a function of their distance to the center. 
This is compared to NFW (a = 1, (3 = 3 and 7 = 1, a = 10 kpc, solid line) and core (a = 0.5, (3 = 3.3, 7 = 0, 
a = 4.5 kpc, dashed line) profiles. 



galactic centre, 7 is the inner slope, (3 is the outer slope, a describes the transition behavior around 
r = a. The popular NFW profile pj] has a = 1, (3 = 3 and 7 = 1 and, adapted to the Milky Way, 
corresponds to , a — 20 kpc and tq = 8 kpc. 

Fig. [2] depicts the radial distribution of the density for all particles in our simulation, calculated by 
the method of Casertano and Hut [33], as will be discussed below. The center is taken as the densest 
point of the galactic halo. Due to resolution limits, this distribution is globally consistent both with a 
NFW-like profile (a = 1, f3 = 3 and 7 = 1, a = 10 kpc) and a cored profile (a — 0.5, (3 — 3.3, 7 = 0, 
a = 4.5 kpc), with ro = R sun = 8 kpc in both cases. At any given radius, the densities exhibit a 
large spread due to non sphericities, local fluctuations and also statistical biases. Moreover, numerous 
density peaks from substructures are also apparent. Notice that the best fit profile is given by the set 
of parameters (a — 0.39, (3 = 3.72, 7 = 0.254, and a = 13.16 kpc), but its physical significance is 
arguable, given the large density fluctuations aforementioned. 

As there is no baryonic component in this simulation, the sun location can be chosen in any direction. 
For the full skymap pictures in the next section, two positions were chosen, one along the positive z 
axis (which is also the projection axis in Figs. Q] and [3]) and one along the positive x axis. The 
corresponding dark matter densities coincide to within 10% and are po = 0.0046 M Q pc -3 = 0.17 
GeVcm -3 and po = 0.0043 M Q pc -3 = 0.165 GeVcm -3 respectively. Throughout the rest of the paper, 
we will use the former for normalization purposes. 

The density around a simulation point i was calculated with the algorithm of Casertano and Hut [33j, 
namely 

P) = ^ P (2) 
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Figure 3: Clump distribution in the galaxy, projected along the z axis. The location of each clump is identified by 
a circle, whose size scales as M 1 / 3 . The red (largest) circle shows the virial radius (as defined in the text) of the 
galaxy (253 kpc). 



where V(rj) = 4-7r/3r 3 is the volume of the smallest sphere around the particle i that includes j 
neighbors. Excluding the particle i itself as well as the jth neighbor in the mass count gives an unbiased 
estimator of the density, with a variance cr| = p 2 /(j — 2). The choice of the number of neighbors used 
to calculate the density is a trade-off between reducing fluctuations and preserving the locality of the 
value computed by this method. A value j ~ 10 seems to give a satisfactory compromise, as we could 
check on some Plummer test models. Higher values of j wash out substructures and inhomogeneities 
present in the data, while smaller values of j imply large statistical uncertainties that mask these 
inhomogeneities with Poissonian noise. As can be seen in Fig. [21 the density inside clumps can be 
several orders of magnitude higher than the density of the smooth component at that location. 



2.2 Clumps 

To identify clumps and subclumps in our simulation, we used the code ADAPTAHOP [34], which 
is an improved algorithm based on HOP that enables to build a tree of structures and substructures. 
Basically, the algorithm divides the simulation points into disconnected groups, or leaves, corresponding 
to local density maxima. To decrease statistical noise, smoothing techniques are applied to calculate the 
density. The connections between leaves are created by performing a search of saddle points between 
groups. The density of a saddle point is then compared with the local maxima on each side, as well as 
a threshold parameter, to decide whether the structures are connected or not. By progressively raising 
this threshold from a minimum value corresponding to a galactic halo overdensity, and performing 
recursively the last check, the algorithm constructs a tree of (sub) structures. We note that the peak 
patches output by ADAPTAHOP are disconnected in space, as they are limited by the closest saddle 
points. As a consequence, clump masses given by this algorithm are often underestimated if other 
structures are present in their neighborhood. We found 108 (sub)clumps attached to the galaxy. Their 
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Figure 4: a) Cumulative subhalo mass function of the clumps in the galaxy obtained with ADAPTAHOP masses 
(red - dark grey) and after the virial mass adjustment procedure (green - lighter grey), b) Cumulative mass for the 
simulation (tot), the smooth component and the clumps, c) A closer look at the cumulative mass profile for the 
clumps (all clumps in red, and only clumps lighter than 10 9 Mq in pink). 



spatial distribution is illustrated in Fig. [31 where each clump is represented by a circle whose radius 
scales as the 1/3 power of its mass, M 1 / 3 . A comparison with Fig. Q] shows that the two pictures are 
consistent. All the clumps within the virial radius are successfully identified by the algorithm. Outside 
this radius, some clumps are found as disconnected from the main galaxy, and therefore not visible on 
Fig. [3 

From this, we derived a cumulative subhalo mass function shown on Fig. 0^. The clumps mass 
fraction inside the virial radius is equal to 3.6% of the total mass within that radius when taking the 
mass values output by ADAPTAHOP. We have corrected these masses by fitting a density profile on 
each clump, and then extrapolating it to the clump virial radius. The adjusted clumps mass fraction 
inside the virial radius is then equal to 5.4%. The mass adjustment has an impact on the mass function, 
as can be seen on Fig. UK- The result is compatible with a power-law with index — 1 above a mass 
threshold of 5 x 10 s M Q . For smaller masses, a flattening of the curve is apparent, due to the resolution 
limit of the simulation, which does not allow the survival of clumps with mass below a given threshold. 
This is around 200 particles, in agreement with [16]. The effect of the mass adjustment further enhances 
this flattening at low masses. Even after this correction, the final mass function that we get is still 
lower than that found by other authors [HI US], but this deviation is not statistically significant. For 
example, we found 3 (7) clumps with mass before (after) adjustment higher than 10 9 M Q) compared 
to 13 in [16]. 

For the clump radial number density distribution, the statistics is too low to obtain reliable informa- 
tion about the behavior near the galactic center (GC). Instead, it is preferable to look at the cumulative 
mass profile as a function of the radius (distance from the GC), shown on Fig. 2b and c. When the 
clumps with a mass > 10 9 M are removed, the cumulative mass profile is well fit by a power-law 
with index n ~ 1.75 up to the virial radius. The corresponding clump mass density profile with index 
n — 3 ~ —1.25 is thus flatter than the smooth component. The mass fraction in the form of clumps is 
therefore increasing with radius. This can be intuitively understood, since clump merging should be 
easier in the central parts. However, the clump density in this simulation is too low to dominate the 
mass fraction at any radius up to the virial radius. 

Examples of density profiles within two clumps are given in Fig. [5] Except for the inner parts, they 
are compatible with a power law with slope ~ —2.5 (FigHJ. It is clear that further in there is a sharp 
transition, but we do not have sufficient points to constrain the inner slopes. In Fig. [5] we have, as an 
illustration, plotted two power laws with slopes and 1, respectively. The transition between the outer 
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slope and the inner slope is described by a concentration parameter c vir = r^ l ir /r_ 2 , where r_ 2 is the 
radius at which d/dr (r 2 p(r)) = and rf ir is the virial radius of the clump. Several models predict 
that the concentration parameter strongly correlates with the virial mass (figure 2 of |23J). 



3 Gamma and neutrino fluxes from dark matter annihilation 

In this work, we consider that the dark matter particle candidate is a typical WIMP provided by some 
new physics beyond the standard model. The gamma ray or neutrino flux per solid angle unit from the 
annihilation of dark matter particles (with mass tudm , density pu m , cross-section (av) , and branching 
ratios BRi into final state i) can be written as 

, M f^ E ^ m l AmM , (3) 

/ v ' 

HEP,, „ =ASTRO 



dQ 4?r 8 m 2 DM Je%-» ^ dE 1>v j^) 



where dNl iV /dE 1<v is the differential gamma/neutrino spectrum per annihilation coming from the 
decay of annihilation products of final state i, the integral is taken along the line of sight with direction 
f2, and 6 = 2 for a self conjugate dark matter particle and 4 otherwise. We have separated into two 
brackets factors that arise from particle physics and from astrophysics. 

The annihilation signal is proportional to the density squared, and can therefore benefit from a 
strong enhancement if the dark matter distribution is highly clumpy. This enhancement is known as 
the boost factor. 



3.1 Astrophysics factor 

To further discuss the enhancement due to the distribution, it is useful to define the dimensionless 
quantity 

j(An) = ^-4r I p 2d -un, (4) 

Ail p^r J 

where the solid angle Af2 can be taken as the experimental solid angle resolution of a given experiment 
such as GLAST. 

To evaluate this quantity for our simulation, two different methods can be used. In the first method, 
p 2 can be calculated at any coordinate of the simulation box space with the Casertano-Hut algorithm. 
Note that the estimator of p 2 is smaller than the square of the estimator for p by a factor (j — 2) /(J — 1). 
The integral along a given line of sight is then calculated with the method of rectangles, with a variable 
step equal to half of the distance to the j = 10th closest neighbor. This ensures that the integral will 
not be overestimated when local density peaks are encountered along the line of sight. Finally, the 
value of J in a cone is the average of the values for different lines of sight within that cone. Fig. [6] 
presents an all-sky view of the astrophysical factor J in a Hammer projection for a value Af2 = 10~ 5 
relevant for GLAST. Since the simulation does not include a baryonic component, the position of the 
observed is only constraied by its distance from the center, i.e. can be anywhere on a sphere of radius 
8 kpc. We calculated the J all sky map for two different viewing positions to allow comparisons. The 
shape of the iso- J contours differs significantly between the two cases and is not circular-like. Indeed 
cosmological simulations show that dark matter halos are not spherical, but have the shape of triaxial 
ellipsoids (e.g. [Ml [36]) and this is true also for the simulation we analyze here. Thus, the observed 
signal will depend on the viewing angle, since the integral in Eq. (J4j) will have different values if the 
integration is e.g. along a major or a minor axis of the ellipsoidal shape. This is also noted in the zoom 
of the central region, shown in Fig. [JJ In this figure we also see that the relevant size of the central 
region is of the order of a kpc, which has been found by galactic dynamic simulations to have a rich 
structure in the baryonic component, including inner bars, inner discs, rings and/or spirals. Whether 
these influence in any way the dark matter in that region still remains to be studied. 
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Figure 5: a) and b) Density profile inside two clumps. The outer slope of the profile corresponds to (3 — 2.5. c) 
Clump luminosity boost factor as a function of the exponent of the clump inner density profile (7) . 



We also evaluated the astrophysical factor with a second method, in which the integral in Eq. (@]) 
is replaced by a finite sum over simulation points. 

The presence of a pole in l~ 2 is potentially dangerous in this method, and can lead to an overestimation 
of the fluxes in case a simulation point is very close to the observer's location. Nevertheless for this 
simulation and for the chosen sun locations, we get comparable results for the two methods, although 
fluxes from low density regions suffer more from statistical noise. To compare the two numerical 
methods more globally, we calculated J and J for the whole sky, Af2 = An. We found a good agreement 
(~ 2%) between the two methods. 



3.2 Particle Physics factor 

The evaluation of the particle physics contribution {HEP ltV ) in Eq. (J3j) is highly dependent on the 
physics beyond the standard model that one assumes. Let us rewrite the HEP term as 

HEP 1>U = I-^Ljv 7i „ , (6) 
Sm 2 DM 



where N 7tV is the number of photons/neutrinos per annihilation 



The spectra coming from dark matter particles annihilation spread up to tudm- The number of 
photons/neutrinos depends on the decay chain of the dark matter annihilation products and these 
annihilations are determined by the particular particle physics model considered. We can estimate 
this number by using an effective and quite model independent approach. Namely, we will assume, 
as in typical BSM model, that dark matter particles (like those coming from SUSY, extra-dimendions 
...) annihilate into SM (Standard Model) particles (W + W~ 7 t + t~ , bb, tt...) whose decays will , in 
turn, produce the gamma/neutrino continuum. The spectrum can then be evaluated from Pythia |37j 
simulations. 
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For gamma rays, the actual number of photons in a given experiment depends on its energy range. 
For GLAST, the energy range will be ~ 1 — 300 GeV. These threshold values fix the actual limits 
for E^ in and E^^. In the 1-300 GeV GLAST energy range, the number of photons estimated from 
Pythia simulations [38] typically yields N 7 ~ 1, 10, 50 and 200 for m DM ~ 10, 100, 1000 GeV and 10 
TeV respectively. 

For neutrinos, we will consider the ANTARES sensitivity and a possible improvement for a km3 size 
telescope located in the mediterranean sea. The typical energy threshold for such a neutrino telescope 
is - 100 GeV leading to 1 < N„ < 10 for 100 GeV < m DM < 10 4 GeV (derived from Pythia neutrino 
spectrum [39]). 

The other important quantity in Eq. J6]) is the annihilation rate (av). From the cosmological point 
of view, (av) can be related to the relic abundance of dark matter, roughly Qdm ~ ^/(av), even 
if strictly speaking the temperature is higher in the primordial plasma at the freeze-out than in the 
surrounding sky at the present epoch and (av)r / (av) freeze-out • For a standard annihilation scenario 
with thermal freeze-out, a value of Qdm like what was found by WMAP implies (av) ~ 10~ 26 cm 3 s _1 . 
Of course, different hypotheses can lead to variations in this value. For instance, one can consider some 
models where coannihilations drive the number of relic particles. Moreover, cosmological scenarios with 
low reheating temperature could allow to decrease the annihilation cross section and still satisfy the 
WMAP constraint. Conversely, scenarios with dark matter particle production through late decays 
(out of equilibrium) of heavier particles could lead to a higher annihilation cross section to respect 
the WMAP relic density. Considering those sources of variability, almost independently of any specific 
particle physics framework, we will assume the reasonable range 10~ 27 < (cru)[cm 3 .s _1 ] < 10~ 24 . 

Consequently, taking into account the production of photons and neutrinos and the annihilation 
rate uncertainties, a typical range for the HEP contribution to the gamma and neutrinos fluxes are 
given by 

10~ 33 < HEP,, [photons cm 3 .s" 1 .GeV" 2 ] < 10~ 26 , (8) 

10~ 34 < HEP V [neutrinos cm 3 .s" 1 .GeV -2 ] < 10~ 28 , (9) 

where the difference between gamma and neutrino is due to the higher energy threshold in neutrino 
experiments, which reduces N v and leads to consider higher values of the dark matter mass0. 

3.3 Comparison with experiments 

With typical values for the local dark matter density po — 0.3 GeVcm -3 [41] and for the Sun to Galactic 
Center distance r$ ~ 8 kpc, one has 

<J> 7 ,„ [photons/neutrinos cm _2 .s _1 ] ~ 10 20 HEP y ^ J AD,. (10) 

3.3.1 GLAST 

GLAST [42] is a satellite which should be launched this year. Its angular resolution should be Ail ~ 
10" 5 srad (i.e. an opening angle of 0.1 degrees ). Considering the sky sensitivity of GLAST given 
by [43], we take 10 -10 photons cm _2 s _1 as a reasonable value to determine the interesting benchmark 
region in excess of which we could expect a signal to be detected by GLAST. 

In Fig. [7j we show a zoom of the central region of Fig. [6] for the two viewing angles analyzed. The 
image size is 60 deg x 60 deg. To smooth the artificial substructures that are due to noise in the 
simulation, we averaged the fluxes in a 5x5=25 pixels square that corresponds to a linear resolution of 
200 pc. Taking into account the HEP 1 contribution, we indicate the J values of 10 3 and 10 2 as quite 
optimistic GLAST benchmarks. The two figures show J values normalized with the same local density. 
One has to keep in mind that the local density can well vary by an order of magnitude depending on 
the Sun location (see Fig. [2]) and that this will influence accordingly the calculated flux value. 

lr Tliis gives typically 2-3 orders of magnitude between EGRET and ANTARES performances despite their similar sensi- 
tivities (~ 10~ 8-7 7(^) cm _2 s _1 ). The same is true for explicit SUSY models [40|. 
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3.3.2 ANTARES and a km3 size neutrino telescope 

AMANDA and ICECUBE are located at the south pole, so that the galactic center region, which is 
the most promising one, is very challenging for those experiments. Thus, we consider in this section 
the ANTARES experiment, whose deployment completion is imminent, and extrapolate the sensitivity 
for a possible km3 size neutrino telescope in the mediterranean sea like the KM3NeT project [44] 
associating ANTARES, NEMO and NESTOR collaborations. The resolution of Antares depends on 
the neutrino-muon angle, but is typically AQ ~ 10 -3 srad. We consider the galactic center sensitivity 
of ANTARES for dark matter derived in [45] , the recent improvement of effective area due to trigger 
update [46] and the expected performance for a km3 size telecope [44]. Taking those references into 
account, we believe that the sensitivity above 100 GeV of a future km3 size neutrino telescope located 
in the mediterranean sea should be around 10 -9 neutrinos cm~ 2 s _1 for point sources in the sky. In 
Fig. [8]we show the central region J skymap but calculated with a Aft — 10~~ 3 resolution corresponding 
to the proposed km3 neutrino telescope resolution. The contour shows an optmistic benchmark region 
corresponding to J — 100 in the zoom for GLAST (no region is available for J = 10 3 ). 

3.3.3 EGRET 

A few years ago, the EGRET collaboration reported an excess in gamma ray fluxes above 1 GeV [47] . 
After subtracting the cosmic ray background, a residual flux of around 5 x 10 _8 photons cm _2 s _1 
remains. The excess is hard to accommodate with natural cosmic ray propagation models, as it requires 
a harder electron injection spectrum or average spectrum in the galaxy different from the local ones [48j . 
As an alternative, the annihilation of dark matter has been suggested. To fit the energy spectrum, the 
dark matter candidate should have a mass around 50 GeV. With an angular resolution Aft — 10~ 3 
srad (i.e. an opening angle of 0.1 degrees), the value of the residual flux would show as the contour on 
Fig. when taking an optimistic HEP 7 that would correspond to J = 10 2 for GLAST. The central 
region would indeed be observable by Egret, while the regions at higher longitudes or latitudes remain 
in the background. However, this is not what was observed by Egret, as the excess is found to be rather 
constant with latitude or longitude, except close to the galactic center region where it is mildly higher. 
When fitted to the data, the hypothesis of dark matter yields a structure which is incompatible with the 
NFW profile. Instead, an prolate isothermal halo supplemented by two rings gives good results [49]. As 
a consequence, a constant boost factor of order 40 is needed in this WIMP candidate scenario, and this 
implies excessive secondary antiprotons fluxes, in clear contradiction with observations [50]. Moreover, 
it has been recently argued that the EGRET excess might be due to a calibration problem rather than 
any unknown astrophysical or exotic contribution |51j . 

3.4 Sensitivity to clumps and resolution 

Taken at face value, our results seem to argue that no clump in the simulation seems to be within the 
reach of GLAST. The intrinsic luminosity of an individual clump is determined by its mass and its 
concentration, but the corresponding flux is damped by the distance square factor. Inside the virial 
radius, the clumps contribute for 8.2% of the total luminosity for a mass fraction of 5.4%, but only 
for 0.2% of the total flux. This, however, depends strongly on the distance from the nearest clump. A 
further point to take into account is the resolution limit of the simulation. Due to that, the innermost 
behavior of the clumps profile in poorly determined, as already mentioned. The intrinsic luminosity of 
a clump can be boosted by increasing the inner slope 7. In Fig. we show the clump luminosity boost 
factor as a function of 7, at constant mass and concentration. For 7 values less than 1.3 or 1.4 the boost 
factor is moderate. In particular, the boost for 7 = 1 (NFW profile) compared to 7 = (core) is only 
around a factor of 2. However, for values of 7 above 1.3 the boost factor increases very spectacularly, 
to reach a a factor of 10 around 7 = 1.4. 

To have a much stronger effect, the concentration parameter should be increased. If we take a NFW 
profile for the clump, it is easy to check that the concentration parameter is directly expressed in terms 
of the scale radius a and the virial radius of the clump r c ; as c V i r = r c ija. Then the total annihilation 
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a) 



b) 



Figure 7: a) Zoom in the 60° x 60° central region of Fig. El i.e. a region centered around the direction of the 
Galactic center. The horizontal bar gives the distance of 1 kpc, which is roughly the region which will be an 
interesting target for GLAST. The contours correspond to J = 10 3 (solid) and J = 10 2 (dotted), b) Same as a) 
but with observer position at (i? sml ,0, 0) instead of (0,0, R sun ) in a). The two plots are normalized with the same 
local density value. 



J resulting from a clump of virial mass Mf ir sitting at a distance I ^> r c \ will be given by 



47r r cl r c i p vir , , . , . 



where p V { r is the virial density (200 times the cosmological matter density), and 



(1 + x) 

for x 3> 1. The clump mass and clump size are related by the following relation 



M cl = A-K Pvir r A cl ln(l + c mr ) - " u (13) 

From this - as p V ir/po — 10 -3 and since the probability of having a clump in our immediate neigh- 
borhood is low - we see that small (in angular size) clumps will be visible by GLAST only if they are 
highly concentrated, c V i r ~ 10 2 — 10 3 . 



4 Conclusion & Perspectives 

We evaluated the gamma and neutrino fluxes from dark matter annihilation in a galactic halo framework 
extracted from a cosmological N-body simulation of the HORIZON project. Although such simulations 
are the most elaborate and realistic framework for this kind of studies, there are very few works 
concerning dark matter detection in N-body frameworks. With reasonable assumptions on the new 
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a) b) 

Figure 8: a) Same as figure[7]with Af2 = 10~ 3 for a km3 neutrino telescope and EGRET, i.e. the 60° x 60° region 
around the Galactic center of the J(Afi = 10~ 3 ) skymap. The light contour corresponds to an hypothetical km3 
size neutrino telescope, the dark one corresponds to the level of the controversial EGRET excess, b) Same as a) 
but with observer position at (R sun ,0, 0) instead of (0,0, R sun ) in a). The two plots are normalized with the same 
local density value. 
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physics beyond the standard model and on the cosmological scenario for dark matter particles, we 
proposed an absolute evaluation of gamma and neutrino fluxes. This allowed to test the galactic 
simulation framework with regard to dark matter indirect detection in future experiments like GLAST 
and the proposed km3 size neutrino telescope in the mediterranean sea extrapolated from ANTARES 
sensitivity and preliminary KM3NeT studies. 

Our framework and even other cosmological simulations have a resolution limit which is lower than 
what would be desirable for this type of studies. Nevertheless, it was possible to reach a number of 
interesting results in our framework. In particular, we showed that the galactic center region is a good 
benchmark part of the sky with regard to the GLAST sensitivity. Even if it is more challenging, this 
region should also be studied by future neutrino telescopes, especially a km3 size telescope. Individual 
clumps stand out clearer in the direction of the galactic anticenter. The concentration parameters of 
the clumps are then the crucial information that will determine the flux. Unfortunately, our simulation 
does not provide enough information on scales below 100 pc to clearly conclude about the detectability 
of individual subhalos. Furthermore, it should be stressed that our results for the central region of 
the Galaxy and for the clumps should be considered as a lower bound. A resolution improvement will 
increase the values of the central density and decrease the size of the central region within which we 
have little or no information. This will increase the signal in the direction of the Galactic center, as well 
as in the directions of the clumps. We can hope that specific geometry, non sphericity and structures 
especially in the central region will be highlighted by future gamma and neutrino observations. 

Cosmological simulations made very important progress in the last few years, and a study of gamma- 
ray induced by dark matter annihilation in the most precise N-BODY framework can be found in [17] 
with a resolution allowing the identification of < 10 6 M sun clumps. Though the quality of this work 
is impressive, this minimal scale is still more than 10 orders of magnitude above the typical WIMP 
free-streaming scale (~ 10~ 12 ~ 4 M sun [52]) and also considerably bigger than the minimal surviving 
clump mass, which anyway is still under debate. 
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